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ROTATION, DIFFUSION, AND OVERSHOOT IN THE SUN: 
EFFECTS ON THE OSCILLATION FREQUENCIES AND THE 

NEUTRINO FLUX 

Brian Chaboyer^, P. Demarque^, D.B. Guenther^ & M.H. Pinsonneault^ 



ABSTRACT 

We have studied the importance of the combined effects of rotation, diffusion, and 
convective overshoot on the p-mode oscillation spectrum and the neutrino flux of the 
standard solar model. To isolate the various physical affects included in the new rotation 
plus diffusion models we also constructed solar models to test the significance of diffusion 
and of overshoot by themselves. In previous studies, models that include helium diffusion 
during solar evolution were found to improve the predicted p-mode frequencies for some 
modes and worsen the agreement for others (Guenther et al. 1993). Here we verify this 
result for both the Bahcall and Loeb (1990) formulation of diffusion and the Proffitt 
and Michaud (1991) formulation of diffusion. We find that the effects of rotation on the 
Sun's structure in the outer layers perturbs the p-mode frequencies only slightly when 
compared to the more substantial effects due to diffusion. In the thin overshoot layer 
(taken here to be 0.1 Hp), we have compared the effect of overmixing in a radiative layer 
versus convective (adiabatic) penetration. Neither radiative overmixing nor adiabatic 
penetration has any significant effect on the |)-modes, probably in part because the 
overshoot layer is constrained to be thin. The predicted neutrino flux in our diffusion 
plus rotation model is 7.12 SNU for CI detectors, 127 SNU for Ga detectors and 5.00 x 
10^ erg cm~^ for the ^B neutrinos; this is approximately half-way between the standard 
solar model without diffusion, and the standard solar model with diffusion alone. 
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1. Introduction 

Standard solar models now predict a p-mode oscil- 
lation spectrum which agrees with the observed oscil- 
lation spectrum of the Sim (Gucnthcr et al. 1992a, b; 
Guzik & Cox 1993), within the estimated uncertain- 
ties of the physical input. This is the result of many 
improvements in the input physics, most significantly 
the advances in astrophysical opacities (Kurucz 1991; 
Rogers & Iglesias 1994). The importance of diffusion 
processes, not normally included in the standard so- 
lar model calculation, have also been explored in so- 
lar models (ProfRtt & Michaud 1991, hereafter PM; 
Gucnthcr et al. 1993; Guzik & Cox 1993; Bahcall 
& Pinsonneault 1992a & b). The diffusion of helium 
from the surface convection zone into the radiative 
layers just below the convection zone was found to 
modify the predicted frequencies of p- modes that are 
most sensitive to the conditions near the base of the 
convection zone (Christcnscn-Dalsgaard et al. 1993; 
Guenther et al. 1993; Guzik & Cox 1993). Specifically 
helium diffusion was found to bring the frequencies of 
p-modcs with ^ = 30 to 50 into closer agreement with 
observations. Furthermore, the depth of the convec- 
tion zone of the solar models that include diffusion is 
closer to the depth derived from inversion of the p- 
mode frequency observations (Christensen-Dalsgaard 
et al. 1991) than standard solar models that do not 
include diffusion. 

The agreement is not perfect, though. Even better 
agreement was found when the efficiency of diffusion 
(Bahcall & Locb, 1990) was reduced by a factor of 
two. In addition, it appears that the improvement 
in the p-modes with ^-values from 30 to 50 is at the 
expense of other p-modcs. That is to say, including 
diffusion affects adversely the agreement for modes 
that penetrate more deeply, ^ = to 30, and modes 
that do not penetrate to the base of the convection 
zone at all, £ = 50 to 100 (Guenther et al. 1993). 

In parallel with our work on diffusion in the Sun, 
we have studied the rotational history of the Sun fol- 
lowing the approach of Endal & Sofia (1978, 1981). 
This work has been carried out using a version of the 
Yale code that includes some of the effects of rotation, 
such as, angular momentum transfer and the associ- 
ated chemical mixing due to rotational instabilities 
(Pinsonneault et al. 1989, 1990). In this approach, 
the radial dependence of the angular velocity is spec- 
ified in convective regions; in the fully convective pre- 
main sequence phase, this plus the total angular mo- 
mentum gives the initial conditions. Angular momen- 
tum is assumed to be lost from the surface convection 
zone of the models due to a magnetic stellar wind. In 



radiative regions, angular momentum is conserved lo- 
cally; the timescales for angular momentum transport 
and mixing are then estimated and solved for, using 
a coupled set of diffusion equations. In such models, 
the spindown of the outer convection zone creates a 
shear layer at its base. 

Rotational mixing will counteract element separa- 
tion; this is especially important below the convection 
zone where both mechanisms arc the most effective. 
In addition, diffusion will interact with other mix- 
ing processes associated with the outward transfer 
of angular momentum from the interior in the Sun. 
This non-linear interaction between rotationally in- 
duced mixing and diffusion has recently been mod- 
eled in a self-consistent way during solar evolution by 
Chaboyer (1993) and Chaboyer et al. (1994) (here- 
after CDP). Earlier work by ProfRtt and Michaud had 
considered the competition between diffusion and a 
hypothetical "turbulent diffusion" law meant to sim- 
ulate rotationally induced turbulence in the Sun. 

The purpose of this paper is to further explore the 
effects of diffusion on the p-mode frequencies. We 
wish to determine if the effort to include a more realis- 
tic treatment of transport phenomena in solar models 
leads to better predictions of the solar p-vaode oscilla- 
tion frequencies. We also wish to investigate the im- 
portance and the nature of the thin transition layer at 
the base of convection zone. In particular, we are in- 
terested in testing whether this transition overshoot 
layer is best represented by simple convective over- 
mixing (in radiative equilibrium), or by convective 
penetration (in adiabatic equilibrium), as discussed 
by Zahn (1991, 1993). 

Finally we estimate the effect of these improve- 
ments on the calculated neutrino flux. Bahcall and 
Pinsonneault (1992a & b) have computed solar mod- 
els with the Yale code to evaluate the expected flux 
of solar neutrinos, and have explored the impact 
of helium diffusion and other changes in the input 
physics for the solar model. In the case of the neu- 
trino flux, the discrepancy between models and ob- 
servations (the "Solar Neutrino Problem") has not 
only remained, but has increased (Turck-Chieze et al. 
1988; Sackmann, Boothroyd and Fowler 1990). Prof- 
fitt (1994) including both helium and heavy element 
diffusion in his calculations, and using a different set 
of diffusion coefficients, has found an even larger in- 
crease in the expected neutrino fluxes. 

The solar models presented here were constructed 
with the rotating version of the Yale stellar evolution 
code to which was added the treatment of thermal dif- 
fusion in the Bahcall and Locb (1990, BL) and Proffitt 
and Michaud (1991) formulations. The constitutive 
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physics differs from the BahcaU-Pinsonneauh (f992a 
& b) models only in some details, such as the treat- 
ment of the Debye-Hii0ckel corrections to the equa- 
tion of state (see Guenther et al. 1992). The models 
arc fully described in Chaboycr (1993) and CDP. In 
section 2, we briefly describe the set of models calcu- 
lated. In section 3 we compare the observed p-mode 
spectrum of the Sun to the standard solar model and 
the altered solar models that include a variety of com- 
binations of the following physical processes: Bahcall 
and Loeb diffusion, Profhtt and Michaud diffusion, 
rotation, and the thin overshoot layer below the con- 
vection zone. The neutrino fluxes are discussed in 
section 4. 

2. The Models 

The Sun is depleted in lithium by a factor of 100 
to 200, and beryllium by 0.3 dex, with respect to 
the meteorites (Anders & Grevcsse 1989). To help 
explain the observed lithium surface abundances, it 
has been suggested that turbulent mixing takes place 
during evolution below the convection zone (Schatz- 
man 1969, 1977; Vauclair et al. 1978; Vauclair 1988). 
Deliyannis (1990) has emphasized the importance of 
explaining both lithium and beryllium depletion si- 
multaneously by the same mixing process. Such mix- 
ing is needed even when the effects of diffusion of light 
elements are taken into account. It has also been 
pointed out that the presence of a fully mixed layer 
at the base of the convection zone, due to convective 
overshoot or penetration, could be detected using so- 
lar seismology (Berthomieu et al. 1993; Monteiro et 
al. 1993; Stix 1994). In rotating solar models (Pin- 
sonneault et al. 1989, 1990), the existence of a mixed 
layer follows naturally from the large rotational shear 
at the base of the convection zone which generates 
turbulence, and in turn mixes the shear layer. In this 
work we examine the significance of diffusion, rota- 
tion, and the mixed layer within the context of the 
p-mode oscillation spectrum. To carry out this inves- 
tigation we have constructed a variety of models that 
include some combination of the effects of rotation, 
diffusion, and a mixed layer. In addition, because 
some authors model the overshoot layer by extending 
the adiabatic temperature gradient below the base 
of the convection zone we also calculate solar mod- 
els with this formulation of overshoot. Specifically 
we have constructed solar models with the following 
characteristics: 

1. ssm — a standard solar model, 

2. nodh — a standard solar model without the 
Debye-Hiickel correction to the equation of state. 



3. adovOl — a standard solar model with a thin 
{Q.lHp) mixed adiabatic layer below the con- 
vection zone, 

4. adov03 — a standard solar model with a thin 
{Q.2>Hp) mixed adiabatic layer below the con- 
vection zone, 

5. mxOl — a standard solar model with a thin 
(0.1 iJp) mixed radiative layer below the con- 
vection zone, 

6. pm — a standard solar model with helium dif- 
fusion in the ProfRtt and Michaud formulation, 

7. pm08 — a standard solar model with helium 
diffusion in the ProfJitt and Michaud formula- 
tion but with the diffusion efficiency reduced by 

20%, 

8. bl — a standard solar model with helium diffu- 
sion in the Bahcall and Loeb formulation, 

9. d-|-r — a standard solar model with rotation 
plus helium diffusion (Proffitt and Michaud for- 
mulation but with the diffusion efficiency re- 
duced by 20%), 

10. d-frZ — a standard solar model with rotation 
plus helium diffusion as in model d-|-r, but with 
the same Z as model pm08. 

The standard solar model and its variations are 
nearly identical to the models described in CDP. Two 
modifications were made to accommodate the p-mode 
frequency calculations: (1) a detailed atmosphere 
based on the Krishna Swamy (1966) empirical fit to 
the solar atmosphere was added to the models, and 
(2) the inner most mass shell was positioned closer to 
the actual center (required for accurate £ = Op-mode 
frequency calculations). 

The models and some of their key distinguish- 
ing characteristics are summarized in Table 1. The 
mixing length and helium abundance of each model 
have been adjusted so that the radii and luminosities 
of the models match each other to one part in 10^ 
(Lq = 3.8515 X 10^3 erg g-i. ^ g 9593 ^ lO^Ocm 
) at the common age of 4.55 Gyr. The solar age from 
the zero-age main sequence has been estimated by 
Guenther (1989, revised in 1992a) to be 4.52 ± 0.03 
Gyr. All the models except the adovOl and adov03 
models were evolved from the fully convective pre- 
main sequence phase. They reached the main se- 
quence in about 0.04 Gyr. 

The models constructed with a mixed adiabatic 
layer extending below the base of the convection zone 
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(adovOl and adovOS) were evolved from a ZAMS 
model (produced during the ssm evolutionary run). 
This introduces a small error of the order of one part 
in 10^ in the temperature and density because of some 
loss of information in the entropy term. The effect on 
the p-mod(^ spectrum is negligible. But it should be 
kept in mind that there is a small systematic differ- 
ence in the calculated neutrino flux of +0.04 ^^Cl SNU 
compared to evolving from the pre-main sequence. 
We did not evolve these models from the pre-main 
sequence because the numerical algorithm to extend 
the adiabatic temperature gradient below the surface 
convection zone lead to convergence problems as the 
model ceased being fully convective. This difficulty is 
only encountered when we impose the high tolerances 
(few parts in 10 '') necessary for solar p-mode calcula- 
tions, and is not present if we relax the tolerances to 
standard stellar values. The question of how to prop- 
erly treat the ad-hoc inclusion of adiabatic overshoot 
during pre-main sequence evolution will be investi- 
gated separately. Here we avoid the question all to- 
gether by starting our evolution for these two models 
on the ZAMS. 

3. The P-Mode Frequencies 
3.1. Introduction 

Diffusion affects the structure of the solar model 
in three ways: (1) it increases the helium abundance 
in the interior and decreases it in the convective en- 
velope; (2) it increases the depth of the convective 
envelope (Guenther 1994) and (3) it changes the sur- 
face value of Z since Z/X is fixed by observation. 
Those p-modes whose inner turning points are near 
the predicted base of the convective envelope, t = 
30 to 50, (see, e.g., Fig.l in Guenther 1994) wiU be 
the most sensitive to the effects that diffusion has on 
the position of the base of the convection zone. The 
frequencies of lower i'-valued p-modes, which pene- 
trate deeper will be affected by the increased interior 
helium abundance, and the higher ^-valued p-modes, 
which probe only inside the convection zone, will be 
affected by the decreased helium abundance of the 
convective envelope. 

The frequencies of p-modcs with ^-values less than 
100 were calculated for the set of models. In Fig- 
ure 1 we show the frequency difference between the 
models and the observations (Libbrecht et al. 1990). 
The mxOl (the standard solar model with an over 
mixing layer extending Q.XHp below the convection 
zone) is not shown because the model is identical in 
all respects to the ssm model. The d-f-rZ model is 
not shown as it was only constructed to test our un- 



derstanding of the effects of rotation and metallicity 
on the neutrino fluxes. Lines connect p-modes with 
common ^-values (not all i?-values are plotted, only a 
representative sample). Bold lines are drawn for £ = 
30, 40, and 50; these modes are very sensitive to the 
conditions near the base of the convection zone. If 
the p-mode spectrum of a model were identical to the 
observed spectrum then the plot would appear as a 
scries of overlapping horizontal lines passing through 

I'model -Vs = 0/xHz. 

The agreement between theory and observation is 
not perfect. Errors in the modeling of the superadia- 
batic layers and the surface boundary conditions are 
responsible for the non-horizontal slope of the lines 
(Guenther 1994). Errors in the surface layers affect 
the frequencies of all modes, primarily by altering the 
average spacing between adjacent in n modes. This 
is seen as a non-zero slope in the frequency difference 
diagram. In the model this error is associated with 
the difficulties of modeling turbulent convection, cal- 
culating low temperature opacities that include the 
effects of molecules, and modeling the surface bound- 
ary condition. For the purposes of this work, we will 
ignore these problems and focus on errors that we can 
actually address. 

Even when we ignore the non-zero slope of the lines 
there still remains the error associated with the scat- 
ter of the lines, i.e., the thickness of the bundle of 
lines. Modes with different ^-values penetrate to dif- 
ferent depths, hence, the fact that lines joining differ- 
ent £- value modes do not lie on top of each other tells 
us that the interior model is not correct to varying de- 
grees at different depths. Improvements to the mod- 
eling of the interior, via opacities, equation of state, 
inclusion of diffusion, etc., should reduce the spread 
or thickness of the bundle of iso-1 p-mode frequency 
lines. 

3.2. Debye-Hiickel Correction and the Signif- 
icance of the Perturbations 

A casual comparison of the thickness of the bundles 
of lines in the plots of Figure J 1 reveal little difference 
among the models. The effects of diffusion, overshoot, 
and rotation are not that large, producing shifts in 
the frequencies of at most ±3 /zHz. The biggest dif- 
ference appears between Figure la and lb. Figure la 
corresponds to our standard solar model, i.e., our ref- 
erence model for this work. Figure lb corresponds to 
a solar model identical to our standard solar model in 
all respects except that the Debye-Hiickel correction 
has not been applied to the equation of state calcu- 
lation. Improvements to the equation of state and 
the opacities have already been shown (Guenther et 
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al. 1992; Christensen-Dalsgaard, D 1990;Christensen- 
Dalsgaard 1988; Ulrich & Rhodes 1983) to lead to im- 
proved p-mode frequency calculations. Here we use 
the model lacking the Debye-Hiickel correction not 
only to emphasize its importance in producing ac- 
curate solar models for p-mode frequency determina- 
tions but also to provide a reference against which we 
will judge the significance of the effects of diffusion, 
overshooting, and rotation on the p-mode frequencies. 

Clearly the effects of diffusion and rotation as in- 
dicated by the rest of the plots in Figure J 1 are not 
as large as the Debyc-Hiickcl improved cqiiation of 
state (the nodh model is the only model that does 
not include the Debye Hiickel correction). Therefore, 
all conclusions that we draw regarding the beneficial 
or detrimental aspects of diffusion, and rotation, as 
implied by the p-mode oscillation spectrum, will not 
be as strongly supported as is the case for the Debye- 
Hiickel correction. 

3.3. Overmixing and Convective Penetration 

We have noted the importance of the existence of 
a thin mixed layer at the base of the convection zone, 
which is needed in the rotation models to describe the 
shear layer. Since there is some qiiestion about the 
effective mean temperature gradient in this layer, we 
have considered two extreme cases. In the first case, 
we simply have allowed the mixed layer to be chemi- 
cally homogenized beyond the base of the convection 
zone by 0.1 Hp, and have let the Schwarzschild cri- 
terion determine whether the local temperature gra- 
dient is radiative or convective (overshooting in the 
terminology of Zahn 1991). Since the opacity is unaf- 
fected by the mixing in this case, the layer is radiative 
(model mxOl). In the second case, we have extended 
the convection zone by applying the convective (adi- 
abatic) temperature gradient by 0.1 Hp and 0.3 Hp 
beyond the formal boundary of the convection zone 
(models adovOl and adovOS, respectively). See Figure 
2. This procedure is often called convective overshoot 
(Ahrens et al. 1992), and is referred to as convective 
penetration by Zahn (1991,1993). 

The model with an artificially mixed layer below 
the base of the convection zone model mxOl) is iden- 
tical, within the numerical accuracy of the model, to 
the standard solar model, and as a result the p-mode 
frequencies are also identical to standard model p- 
mode frequencies. This is expected because the re- 
gion below the solar convection zone is chemically 
identical (except for trace amounts of Li and Be) to 
the convection zone. This is not the case in models 
that include diffusion, where, as a consequence, over- 
mixing can lessen the penetration of He below the 



convection zone. We note that larger overshoot (by 
as much as 0.3 Hp) is required to match the solar Li 
abundance in non-rotating non- diffused solar models 
(Ahrens et al. 1992) by enhanced pre-main sequence 
burning. But overshoot above 0.1 Hp seems ruled on 
other grounds, such as, the observations of Li abun- 
dances in open star clusters which contain young solar 
analogs (Chaboyer et al. 1994b). 

The p-mode frcqiiencies of the solar model in which 
the adiabatic temperature gradient is extended 0.1 Hp 
below the base of the convection zone (adovOl) are 
nearly identical to the p-mode frequencies of the stan- 
dard solar model. Not until the adiabatic overshoot 
reaches 0.3 Hp do we see a significant affect on some 
of the p-modcs. Comparing Figure Id with la we 
see that only the i? = 30 to 50 j3-modes are affected. 
The structure near the base of the convection zone 
(see Figure 2) of the adovOS model is sufficiently dif- 
ferent from the ssm model to perturb the frequen- 
cies of the p-modes that probe this region. We thus 
conclude that we are imable at this time, within the 
present uncertainties, to differentiate between radia- 
tive overmixing and convective penetration of the or- 
der of 0.1 Hp in our models on the basis of the p- 
mode spectrum. Larger amounts of adiabatic over- 
shoot (i.e., > 0.3 Hp) do reveal themselves in the p- 
mode spectrum but are ruled out by the constraints 
set by the observations of lithium destruction in the 
early phases of the evolution of the Sun. Our con- 
clusion differs from that reached by Berthomieu et 
al. (1993), who found an improvement in the p-mode 
spectrum of their convective penetration model over 
the standard solar model. Comparing their "delta fre- 
quency" plots with our plots (Figure 1) reveals that 
the p-modes of their standard reference model are a 
factor of two more discrepant with observation than 
our standard reference model (ssm). We speculate 
that the substantial overshoot layer which they intro- 
duced may have compensated for the shortcomings of 
their reference solar model. This example illustrates 
the pitfalls in interpreting differences between theo- 
retical models and observation when many unknown 
sources of error in the model remain. 

3.4. Diffusion 

Comparing Figures le (model pm). If (model pm08) 
and Ig (model bl) with Figure la (model ssm), it does 
not appear that diffusion improves the solar model 
at all. In fact, the spread of iso-i lines is greater. 
Looking only at the £ = 30, 40, and 50 lines (heavy 
stroked lines) we see that the scatter of these lines is 
reduced. Because the inner turning points of these 
modes are near the base of the convection zone, they 
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are maximally sensitive to the position of the base of 
the convection zone — modes turning below the base 
are maximally perturbed by the abrupt change in the 
temperature gradient while those turning above are 
not. When the frequencies of these modes are com- 
pared to the frequencies of the Sun, an errant location 
in the base of the convection zone in the model is im- 
mediately discerned by the markedly greater spread of 
the p-mode frequency differences. For a model with a 
poorly positioned convection zone base, the iso-£ lines 
for the £ = 30 to 50 p-modes will be tightly bundled 
at low frequencies where all the p-modes have turning 
points above the base of the convection zone and then 
will flay apart {£ = 30 first, then £ = 31, etc.) as the 
frequency increases and the turning points begin to 
dip below the base of the convection zone. 

The location of the base of the convection zone 
in modes that include diffusion of helium, as noted 
by Guenther et al. (1993) and Guzik & Cox (1993), 
more closely agrees with the depth implied by inver- 
sion (Christcnsen-Dalsgaard et al. 1991), hence, so 
do the £ = 30 to 50 p-modes. Unfortunately, not all 
the modes are improved. In fact, the overall scatter 
of the lines, i.e., the bundle thickness for ^ = to 100) 
is increased. There is still something wrong with this 
solar model. Guzik and Cox (1993) have succeeded 
in improving the agreement between model and ob- 
servation beyond that represented here, by adjusting 
the opacities in their models. We believe that errors 
in the opacities just below the convection zone and 
in the outer layers and/or errors in the solar chemi- 
cal abundances (within the present quoted uncertain- 
tics), may be responsible for most of the remaining 
discrepancy. 

The two different formulations of diffusion tested 
here are insignificantly different, with regard to the p- 
mode frequencies. If Figures le, If, and Ig are super- 
imposed on top of each other then slight differences 
can be seen, hut the magnitude of the differences are 
close to the magnitude of the numerical error associ- 
ated with the model and pulsation calculation, which 
is approximately ±0.3/iHz (Guenther et al. 1989). 
The pm08 model is identical to the pm model except 
the diffusion coefficients of the pmOS are scaled to 
80% those in the pm model. We note that the fre- 
quencies of the f = 30 to 50 j3-modes of the pm08 
model are nearly identical to the bl model and that 
the frequencies of the £ = 30 to 50 p-modcs of the 
pm08 and bl models are more tightly bundled than 
the pm model. Because of the numerical uncertain- 
tics in the model calculation we do not attribute any 
significance to this. 



3.5. Rotation Plus Diffusion 

The p-mode frequencies of the models that include 
rotation and diffusion (see Figure Ih; model d-|-r) are 
insignificantly different from the models that include 
just diffusion (compare Figures le, If, and Ig). When 
compared to the standard solar model (Figure la) the 
overall effect on the modes, like that of diffusion by 
itself, is to worsen the agreement. Despite the tremen- 
dous amount of effort required to produce the rotation 
plus diffusion models, we must conclude that rota- 
tion is not the missing ingredient that when included 
in the solar model calculation removes the remaining 
discrepancies between model and observation. Rota- 
tion hardly affects the internal structure and the p- 
mode frequencies at all. 

However, rotation does tend to inhibit the diffu- 
sion, i.e., a model with rotation and diffusion is sim- 
ilar to a pure diffusion model with approximately 
one half the efficiency of diffusion (CDP). By inti- 
mately connecting the effects of rotation with diffu- 
sion (CPD) one perturbs the effects of diffusion on 
the solar model structure. The efficiency of diffusion 
is effectively reduced by the shear induced turbulence 
below the convection zone. Without rotation and the 
associated shear layer, helium diffuses out the base 
of the convection zone and remains within a diffusion 
scale length of the base (0.05 Rq). With rotation, the 
shear induced turbulence below the convection zone 
spreads out the helium that has diffused out of the 
convection zone into deeper layers and it pushes some 
of the helium back into the convection zone, thereby 
reducing the effective efficiency of diffusion. Contrary 
to previous conclusions (Christensen-Dalsgaard et al. 
1993; Guzik & Cox 1993), the presence of a (pos- 
sibly rotationally induced) mixed layer due to tur- 
bulent diffusion just below the convection zone does 
not necessarily affect adversely the calculated p-mode 
spectrum. The effect on the p-modes depends only 
slightly on the adopted formulation of the diffusion 
coefficient. 

4. Rotation Curve 

While the approach of using the p-mode frequen- 
cies described in preceding sections provides a sen- 
sitive test of the effects of internal rotation on so- 
lar structure, it says little about the rotation curve 
itself. This is primarily because rotation is slow in 
these models and barely modifies their internal struc- 
ture. But the distribution of angular momentum 
in the present Sim can be tested directly by obser- 
vations of rotational splittings of oscillation modes 
which penetrate to different depths in the Sun, and 
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with different amplitudes. Recent inversions of the 
lowest 1 modes observed by the IPHIR space mission 
(Toutain & Frohhch 1992), and of the IRIS network 
data (Loudagh et al. 1993) yield a rotation rate at 
0.2 Rq which is higher than the surface value and 
might be as high as 4 times the observed solar surface 
value, which is compatible with some of the models 
constructed by CPD. On the other hand, the inversion 
of intermediate £ p-mode splitting data (Libbrecht et 
al. 1990) yields a flat (solid body) rotation curve 
from the surface to r ~ OARq. Even our minimal 
differential rotation model (see rotation curve in Fig. 
3 and discussion below) predicts a rotation rate in 
the radiative envelope which is too high by at least a 
factor of two when compared to the rotation rate in- 
ferred p-mode splittings in the Sun. This conclusion 
is corroborated by a recent measurement of the solar 
oblateness (Sofia et al. 1994) which puts an indepen- 
dent limit on the amount of angular momentum in the 
outer parts of the solar radiative envelope which is 2 
sigma below the oblateness of the minimal difrerential 
rotation model. 

Another approach which has proved fruitful in test- 
ing models of the evolution of rotating stars is to ap- 
peal to observations of solar analogs at different stages 
of evolution. Observations of the surface rotation pe- 
riods of solar analogs more evolved than the Sun and 
now observed as subgiant stars, are best explained 
if angular momentum buried deep in their interiors 
is being dredged-up into their convcctivc envelopes 
as their convection zones deepen (Pinsonneault et al. 
1989; Demarque & Guenther 1988). There is in ad- 
dition evidence from surface rotation velocities, that 
evolved Sun-like stars have preserved a significant 
fraction of their initial angular momentum in their 
interiors into advanced phases of evolution (Pinson- 
neault et al. 1991). If these objects are true solar 
analogs, we must expect that there is similarly a reser- 
voir of angular momentum in the deep interior of the 
Sun which is now buried in the inner 50% of its radius 
(or 90% of its mass). Some of this angular momentum 
will eventually be dredged up to the surface from the 
deep interior as the Sun evolves into a subgiant. Since 
all other evidence indicates that the outer radiative 
layers of the present Sun rotate slowly, this angular 
momentum must be located within the inner 50% of 
the solar radius. The rotation profiles of the models 
shown in Figure 3 exhibit too much rotation between 
50% and 70% of the radius. 

The minimal differential rotation model shown in 
Fig. 3 was calculated assuming that the inhibitivc ef- 
fects of the mean molecular weight gradient can essen- 
tially be ignored and that the GSF rotational diffusion 



coefficients (James & Kahn 1970,1971) arc ten times 
larger than the default values used in the d+r model. 
The default values give the best fit to the observed 
rotation velocities, and Li abundances of young clus- 
ter stars. This minimal differential rotation model, 
which represents an extreme case, indicates that some 
important physical process or processes, which arc re- 
sponsible for removing some of the angular momen- 
tum from the outer radiative layers of the Sun, are 
missing in our formulation. 

5. Predicted Neutrino Flux 

The structural effects of rotation are small, and so 
rotation by itself does not appreciably alter the pre- 
dicted neutrino fluxes in solar models (cf. Bahcall 
1989). However, recent studies have found that dif- 
fusion of ''He increases the predicted solar neutrino 
fluxes by approximately 11% for the CI and water 
detectors, and 4% for the Ga detectors (Bahcall & 
Pinsonneault 1992a & b; ProfHtt 1994). Rotational 
induced mixing counteracts the effects of diffusion, 
and so models which include the combined effects of 
rotation and diffusion predict neutrino fluxes between 
standard models (no diffusion), and pure diffusion 
models. This is shown in Table 1. We note that the 
neutrino flux of our ssm model is approximately 0.4 
SNU for ^''Cl lower than the best standard solar model 
of Bahcall & Pinsonneault (1992a). Our models use 
the nuclear energy generation of Bahcall (Bahcall & 
Pinsonneault 1992) and the neutrino cross sections 
are also the same as those used in Bahcall & Pinson- 
neault (1992a). The difference in neutrino fluxes in 
the standard models is primarily due to the different 
final ages (4.55 Gyr versus 4.6 Gyr of Bahcall and 
Pinsonneault), which results in a 0.2 SNU decrease, 
and a combination of the different final luminosities 
(3.8515x10^3 crgs-i versus 3.90 x lO^-'' ergs" ^ of Bah- 
call and Pinsonneault) and different helium and heavy 
metal abundances (0.2708 and 0.01880 versus 0.2716 
and 0.01895 of Bahcall and Pinsonneault). There are 
also differences in the treatment of the Debye-Hiickel 
correction (which here includes the effects of electron 
degeneracy), the opacity interpolation procedure, and 
in the chosen values of some of the physical constants. 

The difference in solar neutrino fluxes predicted by 
our standard models and pure diffusion models, 0.8 
SNU increase for ^^Cl, are in good agreement with 
those found by Bahcall & Pinsonneault (1992a & b) 
and Proffitt (1994). The combined models, which in- 
clude rotation and diffusion, predict solar neutrino 
fluxes which are between the standard models and 
the pure diffusion models, with rotation reducing the 
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effect of diffusion on the neutrino fluxes by one-third. 
The last row in Table 1 (model d+rZ) summarizes 
the characteristics of a combined model in which Z 
has been adjusted to be the same as in the pure dif- 
fusion model pm08, with Z = 0.01941. We see that 
0.18 SNU (or approximately half) of the reduction in 
the neutrino flux in d+r is due to the change in Z. 
The rest of the reduction (0.25 SNU) comes from the 
fact that the difi^usion in the core is also being slowed 
down. 

Including the effects of heavy-element diffusion 

does not appreciable effect the predicted p-modc fre- 
quencies, but has a comparable effect to ''He dif- 
fusion on the solar neutrino fluxes (Proffitt 1994). 
Proffitt (1994) determined neutrino fluxes of 9 SNU 
for CI detectors, 137 SNU for Ga detectors and 
6.48 X 10^ cm^ s~^ for the flux for models which in- 
clude ^He and heavy-element diffusion, and suggested 
that turbulent mixing could reduce these amounts by 
a maximum of 0.17 SNU (CI), 0.88 SNU (Ga) and 
0.13 X 10^ cm^ s^^ (^B). Intcr-comparing no diffusion 
models, the pure diffusion model and the combined 
rotation plus diffusion model in Table 1 indicates that 
turbulent mixing due to rotation has a larger effect 
than that suggested by Proffltt (1994). The turbulent 
mixing induced by rotation will lead to reductions of 
0.34 SNU (CI), 1.8 SNU (Ga) and 0.26x10^ cm^ s'^ in 
the predicted solar neutrino fluxes of Proffltt (1994). 

6. Summary 

6.1. Solar internal structure and dynamics 

We have considered the helioseismic properties of 
solar models which take into account rotationally in- 
duced mixing and its interaction with diffusion. We 
emphasize that rotating-diffusive models arc very sim- 
ilar in structure to standard models of the Sun. But 
they differ from the standard models in several subtle 
but important ways (in particular the depth of the 
convection zone and the structure and composition 
profile of the radiative layers just below the convec- 
tion zone) which affects the p-mode frequencies. 

We find that models that include diffusion (either 
the Bahcall-Loeb formulation or the Michaud-Proffitt 
formulation) and models that include diffusion and 
rotation predict p-mode frequencies that are in bet- 
ter agreement with observation for only those modes 
sensitive to the location of the base of the convection 
zone. Higher 1 j3-modes, confined to layers above the 
convection zone, are adversely affected. The effects 
on the p-modc frequencies of rotation cannot be dis- 
tinguished from the effects of diffusion, although the 
latter produces a larger effect. With regard to rota- 



tion, the p-modcs arc only seeing the inhibitive effect 
that rotation has on diffusion via the shear induced 
turbiflent layer, and the diffusion coefficients can be 
adjusted within their error to mimic the effect of rota- 
tion on the p-mode frequencies. At the present time, 
the errors associated with the calculation of the diffu- 
sion coefficients arc too large to enable us to say one 
way or the other whether rotation is a necessary in- 
gredient in the solar model. The diffusion coefficients 
can be adjusted within their error to mimic the effects 
of rotation on the p-mode frequencies. 

On the other hand, as stellar astronomers we can- 
not ignore the considerable success of the rotating 
models in interpreting observations of other Sun-like 
stars of different ages (Pinsonneault et al. 1989, 1990; 
Chaboyer et al. 1994b). These results argue convinc- 
ingly for the existence of a well-mixed transition layer 
in the Sun, most plausibly the result of rotational 
shear. From the y»-mode frequencies, we are unable 
to distinguish between an adiabatic mixed layer (con- 
vectivc penetration) as favored by Zahn (1991,1993), 
and a radiative mixed layer for overshoot layers less 
than 0.1 Fp. 

It is still very difficult to determine the rotation 
proflle deeper inside the Sun. Here the observa- 
tional errors are large, and the models more uncer- 
tain. There is however agreement over the main re- 
sult that the inner part of the Sun rotates faster than 
its surface. This conclusion is in agreement with the 
preliminary inversions from the IPHIR space mission 
(Toutain & Frohlich 1992) and the IRIS network of 
whole-Sun observing stations (Loudagh et al. 1993). 
All observational evidence also points to slow rotation 
in the radiative layers just below the convection zone, 
which indicates that there is a missing angular mo- 
mentum transfer mechanism in our models, which is 
very efficient in fiattening out the solar rotation curve, 
at least in the outer part of the radiative envelope (see 
the rotation profiles in Figure 3). 

Even though the p-modes of our models do argue 
favorably for diffusion plus overshoot, the remaining 
errors in the model (and their affect on the p-mode 
frequencies) are large enough that if they cannot be 
accounted for as errors in the opacity and equation of 
state, then any combination of diffusion, overshoot, 
and overmixing, which is consistent with the observed 
surface abundances, is possible. We do not believe, 
though; this to be the case, and the simplest inter- 
pretation lies in opacity errors. We know that the 
opacities are accurate to only 115% and that smaller 
adjustments to the opacities can correct the p-mode 
spectrum of the solar model. Guzik and Cox (1993) 
have already shown that by tweaking the opacities, 
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the p-mode frequencies of a solar model with diffu- 
sion can be made to agree with the observations to 
almost ±1 yuHz. 

6.2. Solar neutrinos 

Because the rotating-diffusive models only differ in 
small and subtle ways from the standard solar model, 
the predicted solar neutrino flux is not greatly affected 
by the presence of rotation and/or diffusion in the 
Sun. The results are summarized in Table 1. We see 
that the predicted neutrino flux is 6% lower in the ro- 
tation plus diffusion model compared to the diffusion 
only models. 

Rotation by itself does not alter the predicted so- 
lar neutrino fluxes, as the structural effects of rotation 
are small. However, the mixing induced by rotation 
counteracts the effects of diffusion on the solar neu- 
trino fluxes, such that our rotation- diffusion models 
have neutrino fluxes intermediate between the Bah- 
call & Pinsonneault (1992a & b) standard and pure 
^He diffusion models. 
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FIGURE CAPTIONS 



FIG. 1: The frequency differences, model minus ob- 
served (Libbrecht et al. 1990) of a sample of 
p-modes are plotted opposite the observed fre- 
quency of the mode. Lines join common i?- value 
modes for £ = 0, 1, 2, 3, 4, 10, 20, 30, 40, 50, 
60, 70, 80, 90, and 100. The I = 30, 40, and 50 
modes are connected by a thick line. Individual 
plots correspond to the different solar models as 
identified in the text and Tablel. 

FIG. 2: The temperature gradient as a function of 
radius fraction is shown for the ssm, adovOl, 
and adov03 models. 

FIG. 3: The angular rotation velocity is plotted as a 
function of interior mass fraction M^/Mq, Fig 
3(a), and radius fraction r / Rq, Fig 3(b), for the 
specific rotating model used in this paper (d-|-r 
or model VN of Chaboyer 1993), solid line, and 
for a minimal differential rotation model (model 
SF of Chaboyer 1993), dashed line. 
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